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The radial junction (RJ) architecture has proven beneficial for the design of a new generation of high 
performance thin film photovoltaics. We herein carry out a comprehensive modeling of the light 
in-coupling, propagation and absorption profile within RJ thin film cells based on an accurate set of material 
properties extracted from spectroscopic ellipsometry measurements. This has enabled us to understand and 
evaluate the impact of varying several key parameters on the light harvesting in radially formed thin film 
solar cells. We found that the resonance mode absorption and antenna-like light in-coupling behavior in the 
RJ cell cavity can lead to a unique absorption distribution in the absorber that is very different from the 
situation expected in a planar thin film cell, and that has to be taken into account in the design of high 
performance RJ thin film solar cells. When compared to the experimental EQE response of real RJ solar cells, 
this modeling also provides an insightful and powerful tool to resolve the wavelength-dependent 
contributions arising from individual RJ units and/or from strong light trapping due to the presence of the 
RJ cell array. 



Radial junction (RJ) solar cells 1,2 have been fabricated by deposition on top of a matrix of Si nanowires 
(SiNWs), prepared either by template -assisted etching into Si wafers 3 " 5 or by self-assembled growth on low- 
cost substrates 6 " 9 . Prototype devices of this kind indeed demonstrate superior light harvesting performance, 
achieved due to a long optical absorption length while using only a thin absorber/ or junction layer*' 5,10,11 . Though 
several recent works have addressed the full-field simulation of an array of Si nanowire/rods or radial p-n 
junctions 12 " 14 , the light in-coupling and absorption profile in a more sophisticated and complex multilayer p-i- 
n thin film solar cell structure, as well their implications for device performance, remain largely unexplored. 
Within a coaxial RJ thin film solar cell, a sequence of layers of materials with distinctive optical properties are 
deposited around individual SiNW cores 7 ' 91516 . In such a ID cavity-like composite, photonic and resonant mode- 
matching effects in the RJ units via optical antenna coupling or leaky-mode-resonance 13 ' 17 " 19 can play an import- 
ant role in determining light absorption and scattering. These phenomena have been observed experimentally in 
RJ solar cells fabricated over individual upright 20 or surface-lying 15 nanowires. How this resonant absorption or 
scattering will affect the light harvesting behavior in a coaxial p-i-n junction thin film solar cell remains an 
interesting aspect to address. Furthermore, in pursuit of an optimal light harvesting configuration but that can be 
easily electrically contacted, a "forest" of nanowires standing upright on a substrate surface could provide an ideal 
and realistic architecture for making and contacting large areas of RJ solar cells. Therefore, a comprehensive and 
accurate description of the light in- coupling, propagation and absorption profile within the coaxial multilayer 
structure will be very helpful to shed more light upon the critical issues, and to design a new generation of high 
performance RJ thin film solar cells. 

In this study, we focus on a coaxial RJ solar cell, as illustrated in Fig. 1(a), fabricated around SiNWs grown via a 
vapor-liquid- solid (VLS) mechanism catalyzed by a low- melting-point metal such as tin (Sn) 21 " 24 . Low-melting- 
point metal catalysts, including Sn, In, Ga and Bi, allow one to achieve the low-temperature growth (down to 
240°C) of a well-defined matrix of SiNWs 22 ' 25 . This low-temperature process is compatible with becoming a 
mature thin film technology, as well as having the benefit that the metal remnants can be cleaned off easily by a 
simple H 2 plasma etching without the need to break vacuum for ex-situ chemical etching/cleaning 7 ' 9 ' 24 . Figure 1 (b) 
presents a scanning electron microscopy (SEM) view of such a structure. Our recent progress using this structure 
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Figure 1 | (a) Schematic illustrating the coaxial multilayer structure of a radial junction solar cell unit; (b), (c) and (d) show respectively the SEM images 
of a single SiNW, a p-i-n RJ cell and an ITO coated RJ cell. 



has led to the demonstration of high performance radial p-i-n junc- 
tion amorphous silicon (a-Si:H) thin film solar cells with a power 
conversion efficiency of 8.14% 7 , all fabricated in a conventional 
plasma enhanced chemical vapor deposition (PECVD) system. 
Our goal in this work is therefore to establish a comprehensive 
optical model for a coaxially-deployed multilayer RJ thin film solar 
cell, in order to develop critical clues and a new strategy for optim- 
izing such devices. 

Experimental Description 

The RJ solar cells were deposited on a matrix of SiNWs grown on top 
of ZnO:Al coated Corning glass (Cg) substrates. The Sn catalyst 
droplets were formed in situ by a H 2 plasma treatment of a 2 nm 
thick (nominal) Sn layer evaporated on the ZnO:Al substrate. SiNWs 
were grown and doped p-type by flowing a mixture of silane (SiH 4 ) 
and H 2 -diluted TMB [(CH 3 ) 3 B] during the plasma-enhanced VLS 
growth. The detailed fabrication procedure and deposition para- 
meters are available in our previous work 7 . As one can see in 
Fig. 1(b), the typical SiNW features a gradually tapering shape with 
a length of L w — 1.7 |am, a diameter of —40 nm in the middle and a 
sharp tip of <20 nm. As the next step, an intrinsic a-Si:H layer 
(typically —100 nm thick) and an n-type doped a-Si:H emitter layer 
(of — 10 nm) were coated consecutively over the SiNWs, followed by 
a sputtered indium-tin-oxide (ITO) top contact (T ito — 50 nm) over 
the RJ cell. The RJ cells without and with the TCO coating layers are 
shown in Fig. 1(c) and 1(d), respectively, while an internal multilayer 
structure consisting of p-type SiNW-core (p)/intrinsic a-Si:H (/)/n- 
type a-Si:H (n)/TTO contact is depicted in Fig. 1(a). 

Formulation and Simulations 

To carry out an accurate simulation, the optical dispersion n-k curves 
of each material have been extracted from the modeling of spectro- 
scopic ellipsometry measurements on co-deposited planar thin films 
(as shown in S.l in Supplemental Materials). In this RJ simulation, we 
assume the SiNW core is represented by a cylinder of radius R w = 
20 nm and with lengths (L w ) from 0.6 |im to 2.4 |im. The intrinsic a- 
Si:H layer thickness (T,-) has been varied from 20 nm to 140 nm, 
while the thickness of the n-type a-Si:H layer and ITO layers are 
fixed to be typical values of T n = 10 nm and T it0 = 50 nm, respect- 
ively, consistent with what is seen in Fig. 1(b) to 1(d). In our simu- 
lation, a single radial cell is placed at the center of a square simulation 
box (1.5 |im X 1.5 |im) on top of a glass substrate. Polarized light 
propagating along the y-axis originates from the top plane. 
Dissipative scattering conditions are applied to the entirety of the 
wall and bottom boundaries. The finite element solver module in 



COMSOL MULTIPHYSICS was used to carry out the structural 
modeling and optical simulation. 

In contrast to the situation of radial p-n junction solar cells fabri- 
cated on Si wafers 3,26 , the incident light has to go through first the top 
TCO and then the n-type emitter layers before arriving at the inner 
intrinsic a-Si:H absorber in the radial p-i-n junction thin film solar 
cell. As a consequence, assessing the optical propagation and losses 
within the coaxial radial configuration is a critical issue. In the upper 
panel in Fig. 2(a), we show the E y -fie\d distribution in a SiNW RJ 
solar cell unit under different incident wavelengths of X = 350 nm, 
550 nm and 750 nm, with a SiNW length of L w = 1.7 |im and i-layer 
thickness of T { = 100 nm. At X = 350 nm (with high energy 
photons), the incident light is screened by the outer ITO and n-type 
emitter layers, with only a small fraction reaching the buried absor- 
ber i-layer. At long wavelengths (X = 550 nm to 750 nm), the incid- 
ent light propagates deeper into the RJ, and can be better coupled into 
the cavity- like RJ cell, which starts to behave particularly like a wave- 
guide or antenna at X = 750 nm. 

Based on these field distribution results, the local power dissipa- 
tion - that is the local absorption occurring within the RJ multilayer - 
is calculated according to 



PabsW = c-lii gh fcc(l) = c-Ii ight AKk(X)/l 



(1) 



where I Ught is the local electromagnetic energy intensity, c is the speed 
of light in vacuum, a{X) and k(X) are the absorption coefficient and 
the imaginary part of the complex refractive index h(X) = n(X) + 
ik(X). The corresponding power absorption mappings in the RJ cell 
structure are calculated and presented in the lower panel of Fig. 2(b). 
The absorption of high energy photons (2 = 350 nm) mostly takes 
place within the outer ITO and n-type emitter layers, with little 
absorption occuring in the PV active inner intrinsic a-Si:H layer. 
At X = 550 nm, the inner i-layer becomes the dominant zone for 
light harvesting, while the absorption in the c-SiNW core remains 
still very weak, in spite of a strong field within the core region as seen 
in the corresponding E y field distribution in Fig. 2(a). This is because 
at X = 550 nm, the absorption coefficient of a-Si:H is almost two 
orders of magnitude higher than that of c-Si (typically a a . s/;H /a c . s/ — 
10 2 ). At the even longer wavelength of X = 750 nm, when the incid- 
ent photon energy approaches the optical bandgap of the a-Si:H 
absorber matrix, the incident light has been strongly confined within 
the RJ cell. However, the effective absorption realized within the 
intrinsic a-Si:H layer drops dramatically, see Fig. 2(b), while most 
incident photons get absorbed in the center SiNW, which is consid- 
ered as a dead zone because the high doping level in the SiNW core 
electrode. 
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Figure 2 | (a) Upper panels show the simulated electrical field distribution (E y ) through a single RJ cell at different photon wavelengths, while the lower 
panels (b) present their corresponding absorption profile within each material layer in the multilayer RJ structure. 



Close examination of the absorption mapping at X = 550 nm 
shown in Fig. 2(b) reveals yet more intriguing features of the absorp- 
tion profile within a single RJ cell. For example, the highest absorp- 
tion regions in the i-layer are actually located very close to the inner 
i-layer/SiNW interface, as witnessed in the enlarged inset in Fig. 2(b) 
for X = 550 nm, instead of at the outer emitter/i-layer interface, as 
one would expect in a planar p-i-n thin film solar cell. In 
Supplemental Material S2, absorption profiles within the RJ solar 
cell unit at varying incident wavelengths running from 300 nm to 
800 nm are compiled and displayed in an animation clip. We found 
that this feature is a rather common characteristic for wavelengths 
longer than roughly X > 475 nm. Furthermore, to give a quantitative 
estimation of PV performance, we also calculate a solar-spectrum- 
weighted absorption profile P av within a RJ solar cell unit, as shown 
in Fig. 3, defined as 



always preferable to shorten the collection distance of photo -gener- 
ated holes (considering its lower mobility in the a-Si:H matrix). Our 
finding here reveals however a quite different situation, in that a 
higher concentration of photo-carriers are generated in close prox- 
imity to the inner interface and thus a p-type SiNW core should 
become more favorable. This is indeed an important design con- 
straint for optimizing RJ a-Si:H thin film solar cells. 

Furthermore, by integrating the absorption in each material layer 
of the RJ solar cell, we are able to calculate the effective absorption 
power P e f(X) within the absorber i-layer as a function of photon 
wavelength ranging from 300 nm to 800 nm. As we can see in 

Solar-Spectrum-Weighted Absorption Distribution 

I 4.2xl0 4 W/m 3 



Pav =^ Pubs W-Wamu 



(2) 



where W AM i 3 (X) is the weighting factor derived from the standard 
AMI. 5 solar mass spectrum with J Wami.s(^)^ = 1 when integrat- 
ing over X = 300 nm to 800 nm (the most relevant spectrum range 
for a-Si:H absorption). Note that, the x-axis to jy-axis aspect ratio in 
Fig. 3 has been intentionally enlarged to 3 : 1 to deliver a better view of 
the low absorption density variation in the lateral direction. It is clear 
that the absorption intensity profile shows its highest absorption 
zones in the vicinity of the inner i-layer/SiNW core interface, instead 
of the outer emitter/i-layer interface. This phenomenon indicates 
that the incident light has been coupled effectively into the RJ struc- 
ture, following the field distributions of the resonant propagating 
modes within the ID RJ cell cavity. We note that this scenario is 
strikingly different to the situation in planar p-i-n junction a-Si:H 
cells, where the highest absorption zones are always located at the 
outer emitter/i-layer interface, and therefore a p-type emitter layer is 
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Figure 3 | Solar-spectrum-weighted absorption distribution within the 
RJ solar cell unit, over a spectrum through X = 300 nm to 800 nm. 
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Fig. 2(a), the RJ cell can harness the incident light over a cross 
sectional area much larger than its physical diameter, that is, the light 
can be effectively collected into the ID RJ cell (like a waveguide or an 
antenna) 20 ' 27 . This becomes more prominent for the incidence at long 
wavelengths. So, when we hope to address the absorption potential of 
an individual RJ, the finite lateral dimension of the square simulation 
box W sub should be large enough to exclude (or at least minimize) the 
boundary effects. To be sure of this point, we first calculate the P e fj(X) 
curves with different box dimensions, extending from W sub = 
0.8 um to 2.2 um. Their corresponding P e fj(k) curves and overall 
absorption integral (P oa ) from X = 300 nm to 800 nm are presented 
in Fig. 4(a) and 4(b), respectively. We found that the overall integral 
increases gradually with W sub before saturating at W sub =1.5 urn 
and above, indicating that from this point on the RJ cell becomes 
basically optically isolated within the wavelength range investigated 
here. Based on this observation, we define an effective overall absorp 



the thickness up to 140 nm leads to only marginal gain, and only for 
wavelengths >700 nm when the absorption in a-Si:H becomes very 
weak. It is noteworthy that the i-layer thickness of 100 ~ 140 nm 
used here is only 1/3 ~ 1/2 of the typical a-Si:H i-layer thicknesses 
used in optimized planar p-i-n junction a-Si:H solar cell (280 or 
300 nm to balance absorption and collection). A thinner i-layer 
establishes a stronger electric field in the i-layer, facilitating photo - 
carrier separation and collection, and thus helps to improve solar cell 
efficiency and stability 28 ' 31 ' 32 . 

Another important parameter of the RJ thin film solar cell is the 
length L w of SiNW, which represents a new dimension to adjust in 
the 3D RJ architecture, in contrast to their planar counterparts. In 
Fig. 5(c), we plot P ej ^as a function of the length of SiNWs, ranging 
from 0.6 um to 2.4 um. With increasing L w , there is an overall 
enhancement of the i-layer absorption for wavelengths X < 
650 nm, while for the longer wavelengths the absorption enhance- 



tion cross-section area of S oc ~ < fc = (1.5/im) 2 for such a single RJ ment is rather marginal The integral of from 300 nmtoSOO nm 



a-Si:H cell, and extract an enlargement factor (F) of S oc over its 
physical cross section area (S RJ ) of the RJ cell as 



p ^uc 

Srj 



K{R w + Ti + T n + T it0 ) A 2 



(3) 



From this point on, we choose to fix the simulation box dimension to 
W sub = 1.5 um in the following simulations unless otherwise 
specified. 

To understand the optical losses and effective absorption within 
the RJ cell, let us look at the absorption broken down by layer in 
Fig. 5(a), for a typical RJ a-Si:H solar cell unit, as seen in Fig. 1(b)- 
1(d) with L w = 1.7 um, T { = 100 nm and W sub = 1.5 um. The 
absorbed power is plotted the on y-axis in log-scale. It is clear that 
a majority of the absorption losses at short photon wavelengths 
(below 400 nm) are caused by the outer ITO and n-type a-Si:H 
emitter layers. In particular, the n-type a-Si:H accounts for quite a 
large portion of the overall losses, emphasizing the need to adopt a 
wider band-gap emitter, for instance doped a-SiC 28 or uc-SiO x 29 ' 30 
layers. In the long wavelength region, the absorption losses in ITO 
contact become dominant, while the absorption in both i-layer and 
n-type emitter drop rapidly when approaching the bandgap of a- 
Si:H. Figure 4(b) shows the absorption occurring in the i-layer as a 
percentage of the total absorption in the RJ cell, shown for i-layer 
thicknesses (T,-) ranging from 20 nm to 140 nm. As expected, a 
thicker i-layer can boost absorption for longer wavelengths X > 
450 nm, but this tends saturate for T { > 100 nm. Further increasing 



has been calculated and shown in Fig. 5(d), where we see first a fast 
increase of P e jj absorption with the SiNW length, then followed by a 
clear saturation trend. At X > 650 nm, the absorption coefficient of 
a-Si:H matrix drops quickly and the absorption of incident light is 
dominated by the ITO layer and c-SiNW core as shown in Fig. 5(a), 
as well as by the dissipation on the damping boundary walls. To 
quantify the overall light harvesting behavior as a function of the 
length of RJ cell, we fit the data in Fig. 5(d) following the formulation 
of Lambert-Beer law, which can be written as 



i- 



- 0C RJ L w \ 



(4) 



and extract a saturation absorption of P e jp = 2.67 X 10" 12 Wior L w 
— > oo, as shown in Fig. 5(d) by a dashed horizontal line, and an overall 
pseudo-absorption coefficient for such a single RJ cell to be ot RJ = 1.16/ 
/nm (or 1.16 X 10 4 /cm). According to this prediction, a RJ cell with a 
1.7 um long SiNW core already absorbs —86% of the light absorbed 
by an infinitely long RJ cell. This provides a very useful and practical 
guide to evaluate the benefit of growing long SiNWs, against the need 
to keep the SiNW short to facilitate the conformal coating of closely 
spaced SiNWs. 

Furthermore, the impact of the i-layer thickness (T,-) on the total 
light harvesting behavior has been studied. In this simulation, we fix 
the length of the SiNW at L w = 1.7 um. The calculated absorption 
curves for different T { ranging from 20 nm to 140 nm are shown in 
Fig. 6(a). With the increase of i-layer thickness from 20 nm to 80 nm, 
the effective light absorption P e ff enhances greatly over the full 
spectrum, particularly in the range from 450 nm to 650 nm. 
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Figure 4 | (a) Effective absorbed power within the i-layer at specific photon wavelengths and as a function of different simulation box dimension ( W sub ); 
(b) evolution of the integrated absorption within the i-layer. 
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Figure 5 | (a) Absorbed power broken down by material layer, at varying incident photon wavelength for a RJ cell with SiNW length of L w = 1.7 um 
and i-layer thickness of Tj =100 nm. (b) Percentage of effective absorption in the i-layer relative to total absorption in the RJ multilayer structure, 
(c) Evolution of the absorption spectra in the active i-layer with different Si nanowire lengths from L w = 0.6 um to 1.7 (am. (d) full spectrum integrated 
absorption in i-layer as a function of SiNW length. This trend is then fitted by Lambert-Beer law as marked by the red line. 



Interestingly, upon all the P e ff curves, we observe superimposed 
absorption peaks that seem to shift systematically with i-layer thick- 
ness. To clarify this trend, we have normalized all these curves to 
their maxima in Fig. 6(b) to reveal a clear blue- shifting of the absorp- 
tion peaks with decreasing i-layer thickness. To understand this 
point, we have carried out a 2D mode analysis on the cross section 
of a coaxial RJ cell with an i-layer thickness of T { =140 nm. At the 
wavelength of the absorption peak at X = 614 nm, we search for the 
first 6 propagation modes in the RJ cavity and then integrate their 
individual contributions to the power absorbed within the i-layer 
region. We found that the highest contribution comes from a pair 
of propagation modes with mutually orthogonal field distribution 



patterns, one of which is shown in Fig. 6(c). As we can see, for this 
cavity mode, the field distribution is concentrated in the annular i- 
layer, resulting in a maximal useful absorption in the RJ cell. On the 
contrary, other propagation modes found in the RJ cavity (not shown 
here) display effective wave-guiding around the RJ cavity, but have 
their field intensity maxima falling out of the intrinsic absorber layer 
(either in the outer ITO/n-type layer or concentrated into the highly 
doped p-type SiNW core). Such distributions lead to a much lower 
useful absorption for the radial p-i-n junction thin film solar cell. 
This resonant-mode- selection feature/or criterion is indeed a unique 
aspect in the coaxial radial p-i-n cell, as compared to the situation in a 
simple radial p-n junction solar cells. 
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Figure 6 | (a) Absorption spectra in active i-layer while varying i-layer thickness ( Tj) from 20 nm to 140 nm, (b) full spectrum integrated absorption in i 
layer as a function of Tj. (c) field distribution of mode in the RJ cross section with T; =140 nm. 
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Finally, assuming that the absorption of one incident photon (with 
energy higher than the bandgap of a-Si:H) generates only one elec- 
tron-hole pair and a rapid thermalization of the hot photo-carriers 
occurs, the absorption P e ff curve should correspond/or be propor- 
tional to the external quantum efficiency of a single RJ cell (excluding 
recombination issues). We therefore compare, in Fig. 7(a), the nor- 
malized P e ff curves with i-layer thickness of 80 nm (red), 100 nm 
(green) and 120 nm (blue) to the normalized experimental EQE 
spectrum (black), measured on a real RJ cell array shown in the 
top-view SEM image of Fig. 7(b). The i-layer thickness of the RJ cell 
array is around 100 nm with a statistical spread as shown in the inset 
of Fig. 7(c). As a consequence, though there is a specific resonant 
peak position for each i-layer thickness, the weighted- average res- 
ponse curve becomes rather smooth and featureless. As we can see in 
Fig. 7(c), the agreement between the simulated overall absorption 
curve and the actual normalized EQE curve over the short-wave- 
lengths (from X = 300 nm to 550 nm) is indeed excellent, indicating 
that the light harvesting of short-wavelength photons can be mostly 
assigned to the separately contribution from individual RJ cells. In 
other words, the short-wavelength photons have been effectively 
harvested during their first pass through the radial p-i-n junction, 
with basically no contribution from the light trapping among neigh- 
boring RJ cells. At longer wavelengths, we notice a higher EQE res- 
ponse was measured for the actual "forest" of RJ solar cells, compared 
to the simulated response for a single RJ solar cell. This can be seen 
more clearly in Fig. 7(d), where the enhancement ratio of the normal- 
ized experimental EQE over the absorption response of a single RJ 
cell is plotted as a function of incident wavelength. We note that the 
enhancement factor is close to unity for wavelengths below <550 
nm, but scales up monotonically with the increase of the wavelength. 



This observation implies that the contribution of light trapping 
arising from multiple scattering events among individual RJ units is 
indeed wavelength-dependent. This is reasonable if we recall the fact 
that at long wavelengths, the absorption coefficient of a-Si:H absor- 
ber drops rapidly, and therefore allows multiple passes of the incident 
photons through the RJ cell units before being absorbed. Therefore, 
the benefit of strong light trapping gradually manifests itself at longer 
incident wavelengths. Alternatively, the greater experimental long- 
wavelength EQE can also be explained by the overlapping of the 
absorption cross-sections of neighboring RJ cells within a dense RJ 
forest as seen in Fig. 7(b). Though other factors or geometrical con- 
siderations (for example the somewhat random orientation of the 
VLS-grown Si nanowires on glass) could also be responsible for this 
behavior and need to be better understood in future studies, we 
emphasize that this comparison provides a unique perspective, 
emphasizing the relative importance of the light trapping effect in 
a matrix of RJ cells compared with the impact of the geometry of 
individual RJ cells. 

Conclusion 

In summary, we have presented a comprehensive optical model of a 
single radial junction a-Si:H thin film solar cell, which incorporates a 
complete set of accurate optical material properties in each layer. Our 
simulations have been able to reveal the detailed light in-coupling, 
propagation and absorption profile in the complex coaxial multilayer 
solar cell structure, as well as the impact of varying a set of critical 
geometrical parameters. We also investigate the impact of photonic 
effect/resonant-mode enhanced absorption within the quasi- ID cav- 
ity radial junction cells, and NW design criteria for achieving max- 
imal light absorption. These simulation results are further compared 
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Figure 7 | (a) Normalized experimental EQE spectrum (black square-line) measured for RJ cells on VLS-grown SiNWs (example SEM image in (b)), 
compared to the normalized calculated absorption curves for Tj = 80 nm, 100 nm and 120 nm. (c) Normalized experimental EQE and averaged 
simulated absorption curves weighted according to their appearance statistics, (d) Spectrum of the absorption enhancement ratio observed for an array of 
RJ cells (with strong light trapping effect) over that of a single RJ cell, as extracted from (c). 



SCIENTIFIC REPORTS | 4:4357 | DOI: 1 0.1 038/srep04357 



6 



to experimental results on radial junction a-Si:H thin film solar 
cells fabricated on VLS- grown SiNWs. We note that the findings 
developed herein for RJ a-Si:H thin film solar cells are, in principle, 
applicable to radial junction solar cells fabricated by other techniques 
in different material systems, and therefore could provide insightful 
guidelines for structural optimization and improved performance. 

1. Kayes, B. M., Atwater, H. A. & Lewis, N. S. Comparison of the device physics 
principles of planar and radial p-n junction nanorod solar cells. /. Appl. Phys. 97, 
114302 (2005). 

2. Garnett, E. C, Brongersma, M. L., Cui, Y. & McGehee, M. D. Nanowire Solar Cells. 
ANN. REV. MATER. RES. 41, 269-295 (2011). 

3. Lu, Y. & Lai, A. High-Efficiency Ordered Silicon Nano-Conical-Frustum Array 
Solar Cells by Self-Powered Parallel Electron Lithography. Nano Lett. 10, 
4651-4656 (2010). 

4. Garnett, E. & Yang, P. Light Trapping in Silicon Nanowire Solar Cells. Nano Lett. 
10, 1082-1087 (2010). 

5. Sivakov, V. et al. Silicon Nanowire-Based Solar Cells on Glass: Synthesis, Optical 
Properties, and Cell Parameters. Nano Lett. 9, 1549-1554 (2009). 

6. Kendrick, C. E. et al. Radial junction silicon wire array solar cells fabricated by 
gold-catalyzed vapor-liquid-solid growth. Appl. Phys. Lett. 97, 143108-143103 
(2010). 

7. Misra, S., Yu, L., Foldyna, M. & Roca i Cabarrocas, P. High efficiency and stable 
hydrogenated amorphous silicon radial junction solar cells built on VLS -grown 
silicon nanowires. SOL. ENERG. MAT. SOL. C. 118, 90-95 (2013). 

8. Eisenhawer, B., Sill, I. & Falk, F. Radial heteroemitter solar cells based on VLS 
grown silicon nanowires. physica status solidi (a) 210, 695-700 (2013). 

9. Yu, L. et al. Assessing individual radial junction solar cells over millions on VLS- 
grown silicon nanowires. Nanotechnology 24, 275401 (2013). 

10. Peng, K.-Q. & Lee, S.-T. Silicon Nanowires for Photovoltaic Solar Energy 
Conversion. Adv. Mater. 23, 198-215 (2011). 

11. Kelzenberg, M. D. et al. Enhanced absorption and carrier collection in Si wire 
arrays for photovoltaic applications. Nat. Mater 9, 239-244 (2010). 

12. Ningfeng, H., Chenxi, L. & Michelle, L. P. Broadband absorption of 
semiconductor nanowire arrays for photovoltaic applications. Journal of Optics 
14, 024004 (2012). 

13. Kupec, J., Stoop, R. L. & Witzigmann, B. Light absorption and emission in 
nanowire array solar cells. Opt. Express 18, 27589-27605 (2010). 

14. Zhu, J. et al. Optical Absorption Enhancement in Amorphous Silicon Nanowire 
and Nanocone Arrays. Nano Lett. 9, 279-282 (2008). 

15. Kempa, T. J. et al. Coaxial multishell nanowires with high-quality electronic 
interfaces and tunable optical cavities for ultrathin photovoltaics. Proc Natl Acad 
Sci USA 109, 1407-1412 (2012). 

16. Yu, L. et al. Bismuth- Catalyzed and Doped Silicon Nanowires for One- Pump - 
Down Fabrication of Radial Junction Solar Cells. Nano Lett. 12, 4153-4158 
(2012). 

17. Kallel, H. et al. Photoluminescence enhancement of silicon nanocrystals placed in 
the near field of a silicon nanowire. Phys. Rev. B 88, 081302 (2013). 

18. Cao, L. et al. Engineering light absorption in semiconductor nanowire devices. 
Nat. Mater 8, 643-647 (2009). 

19. Cao, L. et al. Excitation of Local Field Enhancement on Silicon Nanowires. Nano 
Lett. 8, 601-605 (2008). 

20. Krogstrup, P. et al. Single-nanowire solar cells beyond the Shockley-Queisser 
limit. Nat Photon 7, 306-310 (2013). 



21. Misra, S., Yu, L., Chen, W. & Roca i Cabarrocas, P. Wetting Layer: The Key Player 
in Plasma-Assisted Silicon Nanowire Growth Mediated by Tin. The Journal of 
Physical Chemistry C 117, 17786-17790 (2013). 

22. Yu, L. et al. Plasma- enhanced low temperature growth of silicon nanowires and 
hierarchical structures by using tin and indium catalysts. Nanotechnology 20, 
225604 (2009). 

23. Yu, L., Alet, P. -J., Picardi, G., Maurin, I. & Roca i Cabarrocas, P. Synthesis, 
morphology and compositional evolution of silicon nanowires directly grown on 
Sn02 substrates. Nanotechnology 19, 485605 (2008). 

24. Rathi, S. J. et al. Tin-Catalyzed Plasma-Assisted Growth of Silicon Nanowires. The 
Journal of Physical Chemistry C 115, 3833-3839 (2011). 

25. Aella, P., Ingole, S., Petuskey, W. T. & Picraux, S. T. Influence of Plasma 
Stimulation on Si Nanowire Nucleation and Orientation Dependence. Adv. 
Mater. 19, 2603-2607 (2007). 

26. Gunawan, O. et al. High performance wire-array silicon solar cells. Prog Photovolt 
Res Appl 19, 307-312 (2011). 

27. Heiss, M. & Fontcuberta i Morral, A. Fundamental limits in the external quantum 
efficiency of single nanowire solar cells. Appl. Phys. Lett. 99, 263102-263103 
(2011). 

28. Green, M. Thin-film solar cells: review of materials, technologies and commercial 
status. J MA TER SCI-MATER EL. 18, 15-19 (2007). 

29. Cuony, P. et al. Silicon Filaments in Silicon Oxide for Next- Generation 
Photovoltaics. Adv. Mater. 24, 1182-1186 (2012). 

30. Fujiwara, H., Kaneko, T. & Kondo, M. Application of hydrogenated amorphous 
silicon oxide layers to c-Si heteroj unction solar cells. Appl. Phys. Lett. 91, 
133508-133503 (2007). 

31. Martin, A. G. Third generation photovoltaics: Ultra-high conversion efficiency at 
low cost. Prog Photovolt Res Appl 9, 123-135 (2001). 

32. Shimizu, T. Staebler-Wronski Effect in Hydrogenated Amorphous Silicon and 
Related Alloy Films. Jpn. J. Appl. Phys. 43, 3257. 



Acknowledgments 

The authors L. YU, J. Wang, S. Qian, J. Xu and Y. Shi thank the following for financial 
support: Jiangsu Province Natural Science Foundation (Young Talent Program No. 
BK20130573) and National Basic Research 973 Program under Grant Nos 2014CB921101, 
2013CB932900 and NSFC under Nos 61036001 and 61204050. 

Author contributions 

L.Y., S.M., J.W. and P.R. designed/conducted the experiments and simulations, wrote the 
main manuscript text and prepared the figures. S.Q., M.F., J.X., Y.S. and E.J. have 
participated in analyzing the data and reviewing this manuscript. 

Additional information 

Supplementary information accompanies this paper at http://www.nature.com/ 
scientificreports 

Competing financial interests: The authors declare no competing financial interests. 

How to cite this article: Yu, L.W. et al. Understanding Light Harvesting in Radial Junction 
Amorphous Silicon Thin Film Solar Cells. Sci. Rep. 4, 4357; DOI:10.1038/srep04357 (2014). 

|>£>k ^ I This work is licensed under a Creative Commons Attribution 3.0 Unported license. 
I^^KS^H To view a copy of this license, visit http://creativecommons.Org/licenses/by/3.0 



SCIENTIFICREPORTS | 4:4357 | DOI: 1 0.1 038/srep04357 



7 



